Lower Bound for the Ground-State Energy Density of a ID Quantum Spin System 
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We present a simple method to calculate systematic lower bounds for the ground-state energy density 
of a ID quantum spin system. 
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The ground- state energy is a quantity of fundamental in- 
terest in the study of quantum spin systems. For example, 
nonanalytic behaviour of the ground-state energy as a pa- 
rameter in the hamiltonian is varied is a canonical signature 
of a quantum phase transition [1]. The ground- state energy 
also occupies a central role in quantum computational com- 
plexity theory: the problem of calculating a good approx- 
imation to the ground-state energy for a class of 2D spin 
systems is complete for the complexity class QMA, which 
is the quantum analogue of NP 12|,|3|,|4J]. 

Obviously, because of the QMA-completeness results of 
Kempe et. al. |3|] and Oliveira and Terhal [4] we expect 
that the ground-state energy is, in general, extremely diffi- 
cult to approximate. However, as these results only pertain 
to highly disordered/frustrated 2D systems, there is some 
hope that it might be possible to efficiently calculate ap- 
proximations to the ground-state energy for ID systems 
and regular 2D systems. 

This hope has been partially vindicated by the develop- 
ment of the density matrix renormalisation group (DMRG) 
(see [0] and references therein for a description of the 
DMRG and relatives). The DMRG provides an apparently 
efficient algorithm to calculate the ground-state energy and 
other local ground-state properties for ID quantum spin 
systems. An exciting extension of the DMRG to 2D quan- 
tum spin systems was recently developed by JfJ. 

Unfortunately the DMRG cannot certify that the estimate 
it provides for the ground-state energy is close to the real 
value. The reason is that there is no way to rule out the pos- 
sibility that the DMRG has become stuck in a local minima. 
For this reason, if a promise that the calculated ground- 
state energy is close to the correct value is required, then it 
is vital to either develop algorithms which provide an esti- 
mate on the distance from optima, or, because the DMRG 
estimates are always upper bounds, provide a method to 
calculate systematic lower bounds for the ground-state en- 
ergy. 

There has been some previous work on lower bounds 
for the ground-state energy: the most important general 
method which has been proposed so far is due to Ander- 
son 10]. This method has been developed further and ap- 
plied to many situations, see eg. JH]. The Anderson bound, 
while nontrivial, is suboptimal; as we'll see later there exist 
systems for which the Anderson bound cannot be system- 
atically improved fast enough to be useful practically. 



In this Letter we describe a simple method to calculate 
systematic lower bounds for the ground-state energy den- 
sity (the ground-state energy per particle) of a ID spin sys- 
tem. For general local ID systems it is an efficient nu- 
merical procedure to extract the lower bound. For some 
systems our method is purely analytic. As a test we apply 
our method to calculate a lower bound for ground-state en- 
ergy density of the XY model. We also apply our method 
to find a lower bound for the average ground-state energy 
density of the disordered heisenberg model. Finally, we 
construct systems for which the Anderson bound cannot 
be improved efficiently, yet our method provides the exact 
answer. 

We consider quantum systems defined on a set of ver- 
tices V with a finite dimensional Hilbert space Ti x , i.e. a 
quantum spin, attached to each vertex x € V. We always 
assume that V is finite owing to the standard difficulties 
in defining a thermal state for infinite quantum spin 
systems. While we take the limit n — > oo we understand 
that this limit is purely formal and, strictly speaking, our 
results pertain only to the situation where n is large but 
finite. 

We will, for the sake of clarity, introduce and describe 
our results for a ring C of n distinguishable spin-| parti- 
cles. Thus, the Hilbert space Tic f° r our system is given by 

■He ®; d y- 

We now introduce the family H n of local hamiltoni- 
ans we are going to focus on. To define our family we'll 
initially fix some two-spin interaction term G which has 
bounded norm: ||G|| < const. (Note that we can, and will, 
accommodate next-nearest neighbour interactions etc. by 
increasing the local dimension of the spins, i.e. by block- 
ing neighbouring spins.) We write the spectral decom- 
position of G as G = X^=o IAjXAjI- By a trivial 
rescaling of the zero point of energy we'll always take 
AminCG) = Ao = 0. Our family H n of local quantum 
systems is then defined by 



Eg 

3=0 



(1) 



where Gj is a translate of G, i.e., it acts nontrivially on 
spins j and j + 1 as G, and as the identity elsewhere. 

Now, it is clear that the ground-state energy eigenvalue 
E > of H will, in general, be strictly positive. In this 
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FIG. 1: Ground-state energy density estimates for H = 
E ie z(l + 7)cJ ® crj+i + (1 - ® <rV +1 versus 7. The blue 
curve is the estimate Eq. d 1 2b applied to the basic even/odd parti- 
tion Eq. 10. The green curve is the estimate Eq. d 1 2b applied to 
the model where blocks of 2 sites have been made. The red curve 
is mean-field theory. The magenta curve is the exact ground-state 
energy density fl4fl . The trivial lower bound 0] that can be de- 
rived by diagonalising the interaction term G in isolation (i.e. by 
pretending each interaction term commutes with every other one) 
is (eo) > -2 for all 7 € [0, 1]. 



case, by using approximate eigenvectors on blocks, it is 
relatively straightforward to argue that E scales as E Q ~ 
c + e n, where c and e are constants. Typically c < 
and e > 0. Because of this scaling it makes sense to talk 
about the ground-state energy density E /n — > eo- It is 
this quantity we wish to bound. 

The central part of our argument relies on the Golden- 
Thompson inequality 111 ill 1211 (this is Corollary IX. 3. 6 in 
H): 



trfe 
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(2) 



where A and B are hermitian operators. To apply the 
Golden-Thompson inequality we divide the hamiltonian H 
into two pieces: 



H = A + B, 



(3) 



where A = 1 G 2j and B = YTj=o ' G 2 j+i. 

Our argument works by first applying the Golden- 
Thompson inequality Eq. © to the thermal state e~^ H us- 
ing the partition Eq. d3j. This gives us the inequality 



n/2-1 1 



Z(J3) = tr(e-^) < tr(e-^e- /3B ) ! 



(4) 



where Z(j3) is the partition function. Next we study the 
expression e~ l3A e~^ B . We begin by writing 



-0G 



M(P) a(3 a a <g> a?, 



(5) 



a, 13=0 



wherea* = [( ?),(? ),(«-),( J ^ )], is the vector 
of Pauli sigma matrices. 



We use this expansion to derive an expression for 



e -PA e -?B 



M((3) aW2 ---M((3) <ia ,o- a o- a ', (6) 

a a = <j a ° <g> cr Ql ® • • • a"— 1 , aj G {0, 1, 2, 3}, < j < 
n — 1, the standard operator basis. 

Substituting the expansion Eq. (|6]l into Eq. (|4| gives us 

tx{e- 0H ) < 2"tr(M(/3) n ). (7) 

The next step is to bound the partition function from below: 

e -pE < tr ( e -/3H) < 2" tr(M(/?)™). (8) 

We now observe, thanks to the positivity of tr(M(/3)™) for 
all n e N, that 



tr(M(/3) n ) < D\\M 



(9) 



where D is the dimension of M. Applying this inequality 
to Eq. d8]l and taking logs gives us 

-(3E <log(D) + nlog(||2M(/?)||). (10) 

After rearranging we obtain the following inequality for the 
ground-state energy density: 

E /n> -log(D)/(/3n)-log(||2M(/3)||)//3, (11) 

for all G [0,oo). 

We will be interested in the large-n limit where n S> 
l/f3. In this case we can ignore the first contribution on the 
RHS of Eq. ([Tit . Thus we obtain our fundamental inequal- 
ity 

E /n>- inf log(||2M(/?)||)//3 (12) 

/3e(0,oo) 

In principle our fundamental inequality Eq. (fT2l l pro- 
vides an analytic lower bound for E /n. In practice, 
however, we need to resort to numerical evaluation of the 
RHS. This is an efficient procedure (in n) and, thanks 
to the continuity of log(||2M(/3)|| 00 )//3 as a function of 
(3 £ (0, 00), provides a certified lower bound for E /n. 

Can the lower bound Eq. ([121 be improved? There are 
at least two ways to proceed. The first method is to com- 
bine I contiguous spins into blocks A and regard each block 
A as a fundamental spin (of dimension 2 l ). We similarly 
block the hamiltonian H, which leads to a new nearest- 
neighbour hamiltonian for the bigger spins. We can then 
apply the procedure we outlined above, albeit using a dif- 
ferent operator basis from standard operator basis. This 
will lead to a new lower bound. It is obvious that this pro- 
cedure cannot lead to a bound which is any worse: after 
all, by taking / = n we recover the exact value of E . 
The second procedure is to look for decompositions of H 

where H = Y2j=o ^3 anc ^ where Hj has a larger mini- 
mum eigenvalue. We then apply the procedure described 
above using Hj instead of Hj. We'll describe a systematic 
procedure to obtain such decompositions in a future paper. 
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There are several obvious generalisations of our method. 
The first generalisation is to calculate e _/3HA for very large 
blocks A of spins via numerical RG methods such as those 
described in fl5ll and lllql . In this case we must put error 
bars on the calculated lower bound because of the trunca- 
tions required by the methods of lfl5ll and iflfill . Apply- 
ing this technique can give us, in principle, arbitrarily good 
lower bounds to the ground- state energy density. 

The second generalisation is to quantum spin systems in 
dimensions higher than one. In this case it is straightfor- 
ward to apply our argument to blocks of spins and thus 
derive an expression for the upper bound which resem- 
bles the tensor contraction patterns investigated in |6|] and 
lfl7ll . Applying the methods described there to approximate 
such tensor contraction patterns will allow us to derive 
lower bounds for the ground state energy of finite higher- 
dimensional spin systems. 

The third generalisation applies to disordered systems. 
It is entirely straightforward to allow the interaction G to 
vary from site to site; the method described above applies 
with essentially no change to such systems. However, in 
this case, we need to average the lower bound over the en- 
semble of possible interactions. In principle this can be 
done analytically for several models. To illustrate the util- 
ity of our lower bound for disordered systems we apply 
this technique to the random antiferromagnetic heisenberg 
model: 

h = y~] Jh&k ■ <Tfc+i, (B) 

fcGZ 

where J k is a random variable with probability distribution 
function /u(x). Applying our bound Eq. (fT2l . appropriately 
modified for disordered systems, and averaging over the 
measure dfi(x) gives us the following lower bound for the 
expected ground- state energy density: 

d»(x) l , 

(14) 

for all (3 e (0,oo), where (A(x)) = J d/i(x)A(x). 

Choosing, for example, dfi(x) = e~ x /TV, where M 
is a normalisation, and choosing x £ [0, oo) provides 
us with the (numerically obtained) lower bound (eo) > 
—0.833610. Compare this lower bound with the lower 
bound obtained from the Anderson bound [7]: (e ) > 
-1.32934. 

Finally, we say a couple of words about the optimality of 
our approach. Consider the interaction term 

G = — 101> (01 1 + 2|10)(10| + 3|11)(11|. (15) 

Now the Anderson bound applied to H = X)j=o yields 
a lower bound for the ground-state energy density given by 
— 1 whereas the method we've developed here yields the 
exact answer: e = 0. Even after blocking m spins the 
Anderson lower bound is still only — 1 /m. 



The decomposition Eq. I© and subsequent derivation 
are strongly reminiscent of DMRG-type methods based 
on matrix product states (MPS) [18]. It is certainly true 
that methods for the ground-state energy based on MPS re- 
duce to variational problems over a restricted class of vec- 
tors in hilbert space and thus must have a well defined la- 
grangian dual which would, at least in principle, provide 
lower bounds for the ground-state energy density. (Further 
investigation of this dual problem will be reported else- 
where). It is currently unclear what, if any, connection our 
method has to this dual problem. 
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